demographics_table <- function(){
  
  table <- read.csv("../release/paper_sumstats/tables/sum_stats/demographics.csv")
  row.names(table) <- c("Single-membered",
                        "Married, no dependent",
                        "Not married, with dependent(s)",
                        "Married, with dependent(s)",
                        "Number of dependents",
                        "HH health severity score",
                        "Income (share of FPL)",
                        "Age",
                        "Over-50",
                        "Living in rating areas 1, 2, or 3")
  
  counts <- read.csv("../release/paper_sumstats/tables/sum_stats/count.csv")
  count_mat <- cbind(counts[, 2], c(NA, NA), c(NA, NA), counts[, 3], c(NA, NA), c(NA, NA))
  row.names(count_mat) <- c("\\quad unique HHs",
                            "\\quad HH-year observations")
  
  MakeLatexTableAdj(
    startcodelist = list(
      "\\begin{table}[H]",
      "\\begin{centering}",
      "\\begin{tabular}{lcccccc}",
      "&  &  &  &  &  & \\tabularnewline",
      "\\hline", 
      "\\hline", 
      "& \\multicolumn{3}{c}{\\textbf{Individual market}} & \\multicolumn{3}{c}{\\textbf{Small group market}}\\tabularnewline",
      "\\cline{2-7} \\cline{3-7} \\cline{4-7} \\cline{5-7} \\cline{6-7} \\cline{7-7}", 
      "\\textbf{Variable} & \\textbf{Mean} & \\textbf{S.D.} & \\textbf{Median} & \\textbf{Mean} & \\textbf{S.D.} & \\textbf{Median}\\tabularnewline",
      "\\hline"), 
    tableinnardslist = list(
      table[, -1], 
      count_mat), 
    tableinnardsroundvec = c(2, 0), 
    inbetweencodelistlist = list(
      list("\\hline", 
           "Number of & & & & & & \\tabularnewline")),
    endcodelist = list(
      "\\hline",
      "\\hline",
      "\\end{tabular}",
      "\\par\\end{centering}",
      "\\caption{Summary statistics on demographic variables} \\label{tab:main_t_sum_stats}",
      "\\end{table}")
  )
}
